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Agile PRT Deconvolution method and systems, and its uses 
Background of the Invention 

1. Field of the invention 

5 This invention relates to radar signal processing. In particular, this 

invention concerns signal processing of agile Pulse Repetition Time (PRT) 
sampled signal transmitted using spread spectrum technique. 

2. Discussion of the background 

10 Radio frequency (RF) pulse-to-pulse agility is a frequency hopping 

known as a spread spectrum technique. By RF pulse-to-pulse agility involves 
changing the carrier frequency per pulse. This implies a much wider 
frequency band than the minimum bandwidth required to transmit and 
receive a pulse. 

15 

The main advantage of RF pulse-to-pulse agility is the increase of the 
resistance against jamming and interference. The spread spectrum has been 
improving radar and communication since the forties and eighties 
respectively. Besides the jamming suppression, such spread spectrum can 
20 also result in a high range resolution, frequency correlation... 

Even when the corresponding Doppler velocity remains constant 
within a burst (i.e. within coherent processing interval), such RF pulse-to- 
pulse agility implies time-varying Doppler spectra. The spectral analysis of 
25 time varying spectra means in radar application Doppler processing when 
Doppler frequency changes from one pulse to another. 

Conventional Doppler processing can not suffice as it is based on 
constant frequency during a burst. This is the reason why the combination of 
30 the RF agility with coherent processing has always been considered difficult if 
not impossible. 

Even when radar signals with the Doppler phase could be extracted 
per pulse, major problems would just begin in Doppler processing, mainly 
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because the received Doppler phase as well as the target radar cross section 
could differ for different carriers. 

In a fully coherent radar, all frequencies are generated from one single 
5 reference oscillator, so the phase coherence is inherent. Carrier frequency 
f RF , intermediate frequencies f A and f L o, sampling frequency f s and pulse 
repetition frequency f PRF are generated from one and only oscillator. In 
addition, RF agile radar requires more than one intermediate frequency f A , n 
so that a number of different frequency per pulse f RF> n can be generated. 

10 

Coherent integration necessitates non-random phases of reflected 
pulses. In a coherent radar, the transmitted phase is known, but such a well- 
controlled phase may be damaged during the propagation. 

15 At time delay t after the nth pulse, the received radar signal s(t n ), 

f n =T+(n+1)fpR T (after analogue-digital conversion), can be modelled as 
follows: 

s(t n ) = a • g(m n - 0) • e mtn) => e Mg 

where a, g() and <K) represent the complex target echo depending on the 
20 target radar cross section, the two-way antenna voltage gain pattern with the 
scan rate co and the target azimuth 0, and the instantaneous phase, 
respectively. The antenna pattern is assumed to be constant within a burst, 
e.g. g()=1. In general, the pulse repetition time t PRT will also be constant 

25 Another drawback of the RF agility is that the RF agility can make the 

target echo a fluctuate independently from pulse to pulse, if a target contains 
many scatterers instead of one single scatterer and, moreover, if the 
individual scatterers move randomly. 

30 If one (dominant) scatterer and no pulse-to-pulse dependence are 

assumed, the phases of reflected pulses remain non-random. Other Doppler 
processing than discrete Fourier transform is needed, unless samples are 
gathered from the pulses with the same frequency. 
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In an extreme target model such as e.g. the Swerling model II, many 
scatterers contribute equally to the echo signal. The resulting amplitude is 
Rayleigh distributed and the resulting phase can only be assumed to be 
uniformly distributed in [0, 2n]. Accordingly, no Doppler processing but 
5 incoherent integration only is applicable. 

By assuming reflectors not interacting or moving randomly, there are 
non-random phases needed. However, the amplitude variability is 
significantly less troublesome than the variability of the phase. In general, it 
10 would add noise in results of the Doppler processing. 

Decorrelation of radar echoes implies no coherence, i.e. no knowledge 
of the phases or, finally, of radial velocities. Fortunately, there are islands of 
frequency correlation that depend on the target dimensions and position but 
15 also on radio frequency and its change. Thus, the RF agility pattern can be 
chosen that result in correlated radar echoes what enable coherent 
processing. 

In any case, time-varying spectra are present in a coherent Doppler 
20 RF agile pulse radar. Thus, the Fourier analysis does not apply anymore. 

In the article "High Accuracy 35 GHz Tracking Radar", Oderland, I., 
Nordlof, Leijon, B., Proceedings of IEEE International Radar Conference 
1990, time-varying spectra were avoided by combining pulses with the same 
25 carrier. This article discloses that such method can be applied only with 
constant pulse repetition time (PRT). However, this waveform lowers Doppler 
range. This means that the range of possible unambiguous Doppler velocities 
decreases fast with increasing number of different carrier frequencies in a 
burst. 

30 

Summary of the Invention 

This invention solves the above-mentioned drawbacks, in particular 
solving the incompatibility between Doppler processing and spread spectrum 
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such providing an improved anti-jamming technique without narrowing the 
Doppler range. 

An object of this invention is a method for deconvolution comprising 
5 the following steps: 

■ [SO] Combining the pulses with the same carrier in a burst, 

■ [S2] Transforming the obtained signals from time to frequency 
domain, 

■ [S6] Deconvolving the obtained spectra 

10 

In a first embodiment of this invention, such a deconvolution method is 
adapted to irregular PRT sampled signal comprising an irregular samples to 
regular zero-padded samples conversion step [S1] between the combination 
and the transformation steps. 

15 

A further object of this invention is a deconvolution system of pulse 
repetition time sampled signal x(t m ) characterised in that it comprises: 

- Means for combining the pulses with the same frequency in a burst, 

- Means for transforming these pulses from time to frequency domain; 
20 - Mean for deconvolving of the spectra. 

in a first embodiment of this further object of the invention, such a 
deconvolution system means for converting irregular samples x{t m ) to regular 
zero-padded samples f(/T f ), these means for converting receives the irregular 
25 pulses grouped by frequency from the means for combining and transmits 
the zero-padded samples to the means for transforming. 

Another object of this invention is an emitting/receiving system using 
such above described deconvolution method comprising an antenna 1, a 
30 reference oscillator 6, an analogue to digital converter 13 and a processor 
which implement the above described deconvolution method. 



Moreover, another object of this invention is the use of such 
deconvolution method in radar system. 

35 
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A further object to this invention is the use of such deconvolution 
method as anti-jamming. 

Brief description of the drawings 

5 

Further features and advantages of the invention will be apparent from 
the following description of examples of embodiments of the invention with 
reference to the drawing, which shows details essential to the invention, and 
from the claims. The individual details may be realised in an embodiment of 
10 the invention either severally or jointly in any combination. 

- Figure 1, block diagram of coherent Doppler RF agile pulse 
radar front-end according to the invention, 

- Figure 2, overlapped bursts with irregular pulse repetition time 
15 samples according to the invention, 

- Figure 3, flow chart of significant steps in the deconvolution 
method according to the invention. 

More Detailed Description 

20 

Figure 1 shows a block-diagram of the transmitter receiver part in a 
radar, which uses the deconvolution method according to the invention. The 
antenna 1of the radar can consist of, for example, a fixed-radiating antenna 
which is connected to a duplexer 2, consisting of, for example, a circulator. 

25 

To the duplexer 2 is connected in a manner known per se a 
transmitting unit 3 and a mixer 10 whose output is connected to the 
intermediate frequency amplifier 11. A local oscillator 6 is connected to the 
mixer 10 via a frequency multiplier 7 and transmits a signal whose frequency 
30 f A , n constitutes an intermediate frequency which is a multiple of the local 
intermediate frequency f L o' fA,n = K^ n ho- 

A pulse modulator 4 is connected to the transmitting unit 3, which 
consist of an RF amplifier, for example, and, indirectly, to the local oscillator 
35 6. The pulse modulator 4 modulates the signal over a certain predetermined 
carrier frequency fRF,n=fA,n+fio with a certain predetermined pulse repetition 
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frequency fp RF =/lo'Kprf. The carrier frequency f RF>n is given by the operator 9 
which adds a first intermediate frequency f LO received directly from the 
reference local oscillator 6 and a second intermediate frequency f Ai „ 
generated by the frequency multiplier 7. The pulse repetition frequency f PRF is 
5 generated by a frequency divider 5 connected to the output of the reference 
local oscillator 6, fp RF =fLo/K PRF . It is made irregular by changing K PRF . 

Furthermore, a signal-treating unit 12 is connected to the output of the 
intermediate frequency amplifier 1 1 and to the reference local oscillator 6. It 

10 consists of a phase detector, which reproduces the Doppler frequency and 
transmits it to an Analogue to Digital (A/D) Converter 13. The Analogue to 
Digital (A/D) Converter 13 is also connected to the output of a frequency 
multiplier 8, which generates a sampling frequency f s . The frequency 
multiplier 8 is connected to the reference local oscillator 6 so the sampling 

15 frequency depend on the intermediate frequency: f s = Ks /Lo- 

The units shown in Figure 1 correspond to those in a coherent Doppler 
RF agile pulse radar front end according to the invention. All frequencies are 
generated from one single reference oscillator, so that the phase coherence 
20 is inherent. The carrier frequency /rf, the intermediate frequencies Ta and /Lo, 
the sampling frequency f s and the pulse repetition frequency f PRF are 
generated from one and only oscillator 6. 

The amplifiers 3 and 1 1 consist of RF amplifier in order that the front 
25 end is an RF radar's one. In addition, the agile radar requires more than one 
intermediate frequency f A>n so than a number of different carrier frequency 
fREn can be generated. 

Figure 2 shows overlapped bursts with irregular pulse repetition time 
30 samples according to a first embodiment of the invention. In a first step of the 
deconvolution method the pulses with same RF are combined as illustrated 
by the hatched RF1 pulses. 
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Irregular sampling has not been used as widely as uniform sampling 
because the time and frequency analyses of irregular samples are rather 
involved. 

5 Random sampling usually involves random numbers added to regular 

sampling times. The DFT-based power spectrum of random samples 
consists of the power spectrum of signal plus additive uncorrected noise. For 
example, even with a noise-free signal sampled at the Nyquist sampling rate, 
the output signal-to-noise ratio could be 1 only. 

10 

Deterministic irregular sampling involves repeating periodically the 
same sequence of K irregular intervals with the mean sampling interval T s . 
The sampling is usually called interlaced when K=2, and multirate or bunch 
when K has an arbitrary length. Since such a sampling set contains K regular 
15 sampling sets, DFT gives K peaks (per one frequency component of a signal) 
within the range limited by the sampling frequency 1/T S . 

For any regular sampling set {nT}, such that T<1/(2B), a real signal x(t) 
limited to a frequency band B, writes as x(t) = J *T dn x(nT).s\nc[2B(t-nT)]. 

20 Irregular samples are well understood in theory, but their method are usually 

too complicated. For an irregular sampling set {t n } t the reconstruction with bi- 
orthogonal bases {sin c[2B(t -nT))} and {¥„(/)}, writes as: 

jc(/) = £ w c tt .sinc[25(/-/ n )] = ]£ ff *(f„ )¥„(/) where c n is the inner product of x(t) 

with {¥„(*)}. If the set {t n } is limited as \t-t n \<\/ (85), ¥„(0 is a Lagrange 

25 interpolation function. Basis functions can also be frames, i.e. bases whose 
orthogonality is not required. The frame conditions are much weaker and 
more useful for practical purpose. 

The NSSL magnitude deconvolution proposed in the American patent 
30 US 6,081,221 is used for the ground radar Doppler processing of interlaced 
sampling scheme. In general, an NSSL sampling set {t m } is multirate with rate 
K and the mean interval Tk- An NSSL time interval (t m+1 - t m ), as well as the 
whole sequence KT S , are integer multiples of the largest common time 
interval T £ , KT S = LT E , so that the smallest regular set {/7V} can contain {t m }, 
35 {t m } c {iT £ }. 
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Thus, the irregular samples x(t m ) are converted to the zero-padded 
samples r(iT £ ) being product of a sampling scheme c,-, c, = 8jT £ - t m ) } and the 
regular samples x{iT £ ). 

5 

Based on this relation: r = diag(o).x (in the vector form), the spectrum 

of x, can be derived as follows: 

dft{r) = #(c) * dft{x) = C • dft(x) 

\dMx)\ = \C\- l .\dft(r)\ 

where C is a Toeplitz matrix whose row vectors are cyclically shifted 
10 dft(c). Since C is singular and, thus not invertible, the NSSL idea is to use the 
magnitudes instead. It is applicable only if there are no complex additions in 
the product Cdft(x) . This condition implies the bandwidth of x, but it is not a 

constraint for most radar. 

15 The spectrum \dft(r)\ contains L replicas of the spectrum \dft(x)\ that 

are weighted by coefficient from |#r(c)| , in a frequency range which is L/K 

times wider than the mean sampling frequency f s . The deconvolution gives 
the strongest replica, i.e. the signal spectrum \dft(x)\ . 

20 Only the NSSL method supports frequencies above the Nyquist 

frequencies, namely up to L/K times the sampling frequency. 

In track radar, where expected Dopplers are reasonably known, the 

signal spectra supported by NSSL may suffice. Namely, using the 

25 magnitudes in 

dft{r) - #(c) * dft{x) = C • dft{x) 

\dMx)\ = \C\'^\dMr)\ 

implies that the spectrum of the signal x cannot be broader than N/L spectral 
lines, where N and L are the number of regular samples (i.e. the length of the 
vectors r and x) and the length of the basic irregular sequence, respectively. 
30 It can also work in the cases where the signal spectrum is broader than N/L, 
but no distance between the spectral lines can be an integer multiple of N/L. 
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The range of unambiguous frequencies may expand limitlessly, but 
clutter filtering becomes more involved with increasing complexity of the 
sampling. 

5 For this purpose, an embodiment of the deconvolution method 

comprises a conversion step. It consists in the irregular samples with the 
same RF to regular zero-padded samples conversion. This is one for each 
radio frequency. 

10 The deconvolution method according to the invention is illustrated by 

Figure 3. This method can be summarised by the following steps: 

[SO: combining step] The pulses are combined by CARRIER (Radio 
frequency) 

[S2: DFT step] The spectrum dft(r) of these regular samples is 
15 computed; 

[S6: deconvolution] The remaining spectra are deconvolved. 

In a first embodiment of the deconvolution method implemented for 
irregular pulse repetition time samples, a further step has been added 
20 between the combining step [SO] and the DFT step [S2]. This step is : 

[S1: conversion step] The irregular samples x(t m ) are converted to 
regular samples r(iT € ), 

The following other steps may exist between the DFT step [S2] and 
25 the deconvolution step [S6]: 

[S3: isolation step] The clutter spectra are isolated by assuming clutter 
spreads over more than a few range gates; 

[S4: estimation] The clutter spectral lines are estimated from the mean 
and the width of the isolated clutter spectra; 
30 [S5: subtraction] The clutter spectra are subtracted from the total 

spectrum dft(r); 

Wherein the steps S5 and S6, the computation are adjusted to the 
clutter type by the clutter spectra given by step S3. 
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In a second embodiment, the spectrum dft(c) can also be computed 
and its L non-zero components found in step S2. In third embodiment, the 
amplitudes of the clutter spectra can be estimated in step S4. So, the 
convolution operations can be reduced, based on the L non-zero spectral 
5 components between the steps S5 and S6. In other embodiments, these 
second and third embodiments can be combined. 



Furthermore, the clutter spectral lines can be estimated from the mean 



in f d =arg[r(r 1 )]/(2OT 1 ) and/or the width in a f = V ln l/>(^)/P( r 2)] <j cl>3af 



10 of the isolated spectra [S4]. Moreover, the amplitudes of the clutter spectra 

r d k Zd]'Z for/>0 
, v • % . n [S4]. 



can be estimated by g (/,*) = 



The convolution operation can be reduced in 

dft(r) = dft(e) * dft{x) = C ■ dft{x) 

\dft(x)\ = \c\- l -\dft(r)\ 

dft(r) = dfi(s) * dft(x) = C - dft(x) 
before the remaining spectra in , , , , are 

^(x^lCl'.j^r)! 

15 deconvolved [S6]. 

One advantage of the deconvolution method according to the 
invention is that it works above the Nyquist frequencies, offers filtering 
method for any kind of clutter, and anti-jamming. 

20 

More generally, such deconvolution system may be used to 
deconvolve any kind of irregular sampled signal using spread spectrum 
technique not only radar one. 



